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Abstract 

Marine invertebrate recruitment may be affected by food limitation during the pelagic larval life stages. In the present study, 
field data on abundance of bivalve larvae along with their prey (small phytoplankton) were examined to see whether they were 
consistent with predictions made by an energetic model of larval requirements. Bivalve larvae were monitored during 2000 at ten 
different study sites in four different areas (Limfjorden, Sylt-Romo bight, Western Wadden Sea and Delta area) along the coast of 
north-western Europe. Calculation of the energetic requirements of the larvae at 15 °C indicated maintenance costs of a 200-pm 
bivalve larva to be 1.9 x 10“^ J larva” * d” *, while the maximum assimilation rate, resulting in maximum growth, would amount 
to 6.2 X 10”^ J larva” * d” *. Calculation of potential assimilation rates of larvae in the field resulted in estimates between 10”^ and 
10”^ J larva”'d” \ Maximum larval concentrations in the field occurred from May to September and ranged between 17 and 392 
larvae dm”^. Most larvae were able to cover their maintenance costs, but not to attain maximum growth rates. Between April and 
September, the potential assimilation rate averaged 7-26% of the maximum assimilation rate. Under the assumptions made for the 
present study, it is suggested that growth of larvae in north-west European waters is often food-limited. 

© 2005 Elsevier B.V. All rights reserved. 
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1. Introduction 

For most soft-bottom invertebrate species, it is not 
clear which larval processes are crucial for recruitment 
success (reviews in Rumrill, 1990; Olafsson et al., 
1994; Menge, 2000). Among the first to recognise 
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that food limitation of larvae may affect marine inver¬ 
tebrate recruitment was Thorson (1950). In his view, 
larvae could die directly from starvation, or indirectly 
from predators due to a prolonged planktonic develop¬ 
ment time. Experimental evidence has demonstrated 
that food limitation of invertebrate larvae is more likely 
at higher trophic levels (crustaceans) than at lower 
trophic levels (molluscs, echinoids) (Hairston et al., 
1960; review in Olson and Olson, 1989). Frufhermore, 
food limitation is expected to be more important in 
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oceanic waters than in coastal waters, due to the lower 
primary production (Huntley and Boyd, 1984) that 
would result in a lower per capita larval food intake 
rate (Holling, 1959). However, larval food limitation 
has also been reported in coastal waters (Fenaux et al., 
1994; Hansen, 1999). 

To assess whether food limitation is common to 
bivalve larvae, it is necessary to study both larval 
concentrations and concentrations of their prey, prefer¬ 
ably in time and space (Olson and Olson, 1989). Time 
series of larval concentrations in north-western Europe 
have been provided by several authors (Pulfrich, 1997; 
Gunther et al., 1998; De Vooys, 1999; Strasser and 
Gunther, 2001). However, no extensive geographical 
comparison of larval concentrations has yet been made. 
Assuming that the areas may be regarded as closed 
systems and that the timing of spawning of several 
bivalve species in north-western Europe depends on 
species-specific threshold water temperatures (Mytilus 
eduUs: Wilson and Seed, 1974; Podniesinki and McA- 
lice, 1986; Thorarinsdottir, 1986; Crassostrea gigas: 
Ruiz et al., 1992; Macoma balthica: Lammens, 1967; 
Drent and Luttikhuizen, ms. in prep.), it is expected that 
in a comparison, the differences in the height, duration 
and frequency of larval peaks (all bivalve species in¬ 
cluded), between areas primarily depend on (1) the 
spatial uniformity/irregularity of water temperature in¬ 
crease, (2) the number of species, (3) the adult stock 
size, and (4) total water volume. 

Prey concentrations for bivalve larvae have not ex¬ 
plicitly been determined, due to the limited information 
on natural diets of larvae (Baldwin and Newell, 1991, 
1995). In general, the major food source for marine 
invertebrate larvae is considered to be phytoplankton 
(review in Strathmann, 1987), but other sources such as 
small heterotrophic flagellates, amino acids (Langdon, 
1983; Manahan, 1983, 1989), bacteria (Baldwin and 
Newell, 1991) and dissolved organic matter (DOM), 
can also be ingested by larvae. Utilisation of DOM 
(osmotrophy), however, probably only supplements 
nutrients obtained by other feeding modes (Olson and 
Olson, 1989). 

To relate food availability to growth of organisms, 
the Dynamic Energy Budget model (DEB-model; 
Kooijman, 2000) can be used. The DEB-model quan¬ 
titatively describes energy flows through organisms and 
the changes in these flows in environments in which 
food concentrations and temperatures vary. In this 
model, all physiological and ecological processes are 
related in a systematic and coherent way. In short, 
organisms ingest energy, which is assimilated and di¬ 
rected towards growth and maintenance of structmal 


cells, and towards the production of cells used for 
maturation/reproduction and their maintenance. 

The aims of the present study were (1) to compare 
abundances of pelagic bivalve larvae along the coast of 
north-western Europe (Lim^orden, Sylt-Romo bight, 
western Wadden Sea, Delta area) in one year (2000); 

(2) to estimate to which degree bivalve larvae were 
food-limited at these study sites. 

2. Materials and methods 

2.1. Study sites and sampling methodology 

Bivalve larvae were sampled in 2000 in four areas 
along the coast of north-western Europe: (1) Lim^or- 
den in the north of Denmark, (2) the Sylt-Romo bight, 

(3) the western Wadden Sea and (4) the Delta area (Fig. 
1). Samples were collected routinely from jetties or 
piers at high tide or close to high tide, and from boats 
at fixed times irrespective of the tide. Bivalve larvae 
sampling frequency was 20-252 y“* (Table 1). 

Larvae in different areas were sampled with compa¬ 
rable methods (Table 1). Generally, a sample of 10-100 
dm^ of surface seawater was taken with a bucket or 
pumped and filtered over a mesh size of 40-80 pm. 
Samples were preserved in formaldehyde or lugol and 
stored in the dark until further analysis. Only the total 
numbers of larvae were coimted, irrespective of species 
composition, since microscopical discrimination be¬ 
tween larval species is complicated (Loosanoff et al., 
1966; Chanley and Andrews, 1971; Le Pennec, 1980; 
Pulfrich, 1997) and discrimination on the basis of DNA 
analysis in our study region is still under development 
(Van Bleijswijk, ms. in prep.). However, based on the 
stock biomass and on species determination of bivalve 
larvae in earlier years, the majority of the larvae prob¬ 
ably consisted of Mytilus edulis, Cerastoderma edule, 
Ensis americanus, Crassostrea gigas, Macoma 
balthica and Mya arenaria (De Vooys, 1999; Strasser 
and Gunther, 2001; Dolmer and Frandsen, 2002; Smaal 
et al., 2002; Ysebaert et al., 2003). 

Phytoplankton data for the Lim^orden stations were 
collected simultaneously with the bivalve larvae. Phy¬ 
toplankton data for the other stations were obtained 
from phytoplankton monitoring programmes by J. van 
Beusekom (AWI-Sylt), G.C. Cadee (Royal NIOZ), and 
the RIKZ (Middelburg, the Netherlands). Sampling in 
the Oosterschelde was at locations close to larval 
sampling locations (0-5 km). Phytoplankton data for 
stations OS-C and WS were not available. Therefore, 
only data on larval abimdance are shown for these 
stations. 
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Limfjorden 



Fig, 1. Location and abbreviations of names of sampling stations (#) in four coastal waters of north-western Europe. 


2.2. Larval food availability 

Phytoplankton size estimations were based on Dodge 
(1982), Round et al. (1990), Tomas (1997), Kuylen- 
stiema and Karlson (1998), Homer (2000), and Ichise 
et al. (2004). Only phytoplankton taxa <20 pm were 
included in further analyses, based on the assumptions 


(1) that larvae only ingest food particles up to the diam¬ 
eter of their mouth (20.5 ± 1.5 pm SD for M. balthica 
larvae of 184.9 ±9.2 pm, n=15; Van Iperen, pers. 
comm., 2003), which is confirmed by nutritional studies 
(Riisgard and Randlov, 1980; Hurley et al, 1997; Bald¬ 
win and Newell, 1991, 1995) and (2) that the average 
size of larvae in the field is -200 pm (Pulfrich, 1997). 


Table 1 

Overview of methods used to sample bivalve larvae at sampling stations along the coast of NW Europe in 2000 


Area 

Station 

Sampling 

Sampling 

Interval 

Station 

Moment 

of 

sampling 

Sampling 

method 

Mesh 

(pm) 

Sample 

Volume 

(1) 

Fixation 

Method 

Analysed 
Volume 
(% of total) 

Limfjorden 

LF-LB 

Jan-Dec 

1^ mo^ ‘ 

Boat 

- 

PP 

60 

450-1500 

L 

0.2-100 


LF-SF 

Jan-Dec 

l^mo-‘ 

Boat 

- 

PP 

60 

450-1500 

L 

0.2-100 


LF-NB 

Jan-Dec 

1^ mo^ ‘ 

Boat 

- 

PP 

60 

450-1500 

L 

0.2-100 

Sylt-Romo bight 

Sylt 

Jan-Dec 

1-7 w^ ‘ 

Ferry pier 

HT-1 h 

Bu 

80 

10 

F, H 

100 

Western 

WWS 

Jan-Dec 

1-3 w^‘ 

Jetty 

HT 

Bu 

40 

10 

F, B, R 

5 

Wadden Sea 
Delta area 

OS-N 

Apr-Oct 

1 w^‘ 

Boat 


Bu 

55 

100 

F, B 

13 


OS-E 

Apr-Oct 

1 w^‘ 

Boat 

- 

Bu 

55 

100 

F, B 

13 


OS-C 

Apr-Oct 

1 w^‘ 

Boat 

- 

Bu 

55 

100 

F, B 

13 


OS-W 

Apr-Oct 

1 w^‘ 

Boat 

- 

Bu 

55 

100 

F, B 

13 


ws 

Apr-Oct 

Weekly (moon) 

Jetty 

HT-1/2 h 

Bu 

55 

100 

F, H 

10 


HT=high tide, LT=low tide, Bu=sampling with 10-litre bucket at surface, PP=sampling with plankton pump (150 1 min^') at 1 m depth, 
F=preserved with 4-10% formaline, H=sample buffered with hexamethylenetetramin, B=buffered with borax, R=coloured with Bengal Red, 
L=preserved with Lugol. 
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Phytoplankton taxa forming chains, colonies or posses¬ 
sing large spines were also excluded from the analysis. 

Biovolumes of the selected phytoplankton species 
were obtained from the literature or estimated using a 
standardised set of equations for biovolume calculations 
(Hillebrand et al., 1999). Biovolumes were converted to 
carbon content (C; pgC pm“^), using the conversion 
logioC = -0.422-l-0.758 logio V for diatoms (Strath- 
mann, 1967), logioC = -0.363 -I- 0.863 logio V for 
other nanoplankton, including dinoflagellates (2-20 
pm; Verity et al., 1992) and C = 0.165 V for ciliates 
(Putt and Stoecker, 1989), where V (pm^) is the volume 
of the cell. Subsequently, carbon content was converted 
into energy content (E; J cm“ ^) using the conversion 1 g 
C=46 X 10^ J (Salonen et al., 1976). Next, cell carbon 
content was multiplied by cell concentration and 
summed for all species per sampling date. 

2.3. Dynamic Energy Budget model 



Fig. 2. Energy flow through an organism in the DEB-model, after Van 
Haren and Kooijman (1993). Rates: 1 ingestion, 2 defecation, 3 
assimilation, 4 mobilisation of energy into reserves, 5 demobilisation 
of energy fi-om reserves, 6 utilisation, 7 growth, 8 maintenanee, 9 
maturation maintenance, 10 maturation. The rounded boxes indicate 
sources or sinks, the squares indicate state variables. 


Calculations of maximum assimilation rate (J 
larva” * d“'), assimilation rates as a function of food 
concentration (J larva” ^ d” ^) and maintenance costs (J 
larva”'d”') were based on predictions of the Dynamic 
Energy Budget (DEB) model developed by Kooijman 
(2000) and adapted for M balthica larvae by Philippart 
et al. (ms. in prep.) (Table 2). The DEB-model 
describes the energy flows through an organism (Fig. 
2) and the changes in these flows in environments in 
which food concentrations and temperatures vary. 


In the model, energy intake Jx is described by a 
type-II functional response (Holling, 1959) (Fig. 3), in 
which the ingestion rate depends hyperbolically on the 
available food concentration: 


Jx = {Jxm} ■ 


X 




where Jx is the ingestion rate (cells larva”* d”*; for 
other symbols see Table 2). The ingested food is con¬ 
verted into energy reserves with constant assimilation 


Table 2 

Description of symbols in Dynamic Energy Budget (DEB) model and parameter estimates and eonversion faetors for the DEB model of a Macoma 
balthica larva at 15°C (after Philippart et al., ms. in prep.) 


Factor 

Description 

Value 

Ref 

Units 

X 

Food density 

- 


cells cm“^or J cm“^ 

Xk 

Half-saturation constant for feeding 

10,000 

1 

cells Isochrysis cm“ ^ 

{Jxm} 

Surfaee-specific maximum ingestion rate 

3.06 


cells Isochrysis 

{Pam} 

Surfaee-specific maximum assimilation rate 

1.04 


J 

6m 

Shape coeffieient 

0.386 


- 

K 

Energy allocation coefficient 

0.76 


- 

[Ed 

[Pm] 

Volume-specific costs of growth 

Volume-specific costs of somatic maintenance 

5.67 X 10^* 

3.30 X 10-" 

2 


[Pj] 

Volume-specific costs of maturity maintenance 

3.30 X 10-" 


J |xni“^ 

V 

Volume at 200 pm shell length 

(0.386x200)^ 


\im^ 

Ta 

Arrhenius temperature 

7596 

3 

K 

k(Ti) 

Physiological rate at reference temperature Ti 

- 


- 

Ti 

Reference temperature 

288 (15) 


KCO 

k(T) 

Physiological rate at temperature T 

- 


- 

T 

Seawater temperature 

” 


K 

Conversions 


Energy content of Isochrysis cells 

6.07 X 10-’ 

4 

J cell ' 


Assimilation efficiency coefficient 

0.56 

3 

- 


(1) Jespersen and Olsen, 1982; (2) Beukema and De Bmin, 1979; (3) Van Haren and Kooijman, 1993; (4) Sprung, 1983. 
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j* 0.005 

e 

^ 0.000 

I 
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0.010 
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Fig. 3. DEB-model predictions of energy intake (J larva^' d^*) as a function of food concentration (functional response; cells Isochrysis 
galbana cm“^) of a 200 pm Macoma balthica larva at 4 different temperatures (0, 10, 15 and 20°C; Philippart et al., ms. in prep). From top 
to bottom, lines indicate maximum ingestion rate (Jxm). ingestion rate (Jx)> maximum assimilation rate (pAm). assimilation rate (pa), and 
maintenance costs (M). The vertical line indicates the half-saturation constant X^. 




0 20 40 60 80 100 0 20 40 60 80 100 

Food concentration X (10^ cells Isochrysis ml'’) 


efficiency. The assimilation rate can therefore be de¬ 
scribed by: 


PA = {pAm} ■ 


- ■ 


where Pa= assimilation rate (J larva”* d“*; for other 
symbols see Table 2). 

The assimilated energy flows into a storage com¬ 
partment (Fig. 2) and no maintenance costs are paid 
over storage of energy. From this compartment, the 
fraction k of the mobilised reserve is mobilised into 
somatic growth and maintenance of somatic tissue, 
while the remainder (1 -k) is converted into growth of 
maturation tissue and maintenance of this tissue. 

Maintenance costs of both somatic and maturation 
tissue are proportional to volume, since each cell uses a 
fixed amount of energy per unit of time. Somatic 
maintenance has priority over growth, hence an organ¬ 
ism stops growing as soon as the fraction k of the 
mobilised reserve is equal to, or lower than, somatic 
maintenance. In this model, somatic maintenance (Fig. 
2) is expressed as: 


PM = \pm] ■ V, 


where Pm= somatic maintenance costs (J larva * d *; 
for other symbols see Table 2). 


Maturation maintenance (Fig. 2) is expressed as: 

where Pj=maturity maintenance costs (J larva”* d”*; 
for other symbols see Table 2). Together, maturity 
maintenance costs pj and somatic maintenance costs 
Pm are referred to as maintenance costs M (J larva”* 
d”*) in the present study. 

Temperature affects all physiological rates (Kooij- 
man, 2000; Fig. 3). In DEB, the Arrhenius description 
has been selected to describe the temperature correction 
for a physiological rate k(T) at a particular temperature 
T according to: 


k{T) = kTi ■ 


where k(T) is the physiological rate at temperature T(K; 
for other symbols see Table 2). 

In the present, study the length of the M balthica 
larva in the DEB-model (Philippart et al., ms. in prep.) 
was set at 200 pm, as an approximate average length 
of larvae found in the water column (Pulfrich, 1997). It 
is assumed that the energy budget of other bivalve 
larvae of the same length is in the same order of 
magnitude. 
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In the DEB-model for M balthica larvae (Philippart 
et al., ms. in prep.), the half-saturation constant was 
based on the ingestion rates of mussel larva at 15 °C 
(Mytilus edulis\ Jespersen and Olsen, 1982) and 
set at 10000 Isochrysis galbana cells cm“^ (Table 2). 
Assuming an energy content of I. galbana of 6.07 x 
10“^ J cell”' (Sprung, 1983), and an area-specific 
maximum ingestion rate {Jxm} of 3.06 I. galbana 
cells pm“^ d“' (Jespersen and Olsen, 1982), the area- 
specific maximum ingestion rate {Jxm} could be ex¬ 
pressed as 3.06 X 6.07 X 10”’=2.5 x 10“® Jpm”^ d“'. 
The assimilation efficiency was assumed to be 0.56 (Van 
Haren and Kooijman, 1993), yielding an area-specific 
maximum assimilation rate {pAm} of 2.5xl0”®x 
0.56=1.4x10“® J pm“^ d“'. Somatic maintenance 
costs [pm] were calculated from the decrease of somatic 
mass during the non-feeding period of adult M. balthica 
(Beukema and De Bruin, 1979), and were 3.30 x 10“ " 


J pm“^. The maturity maintenance costs [pj] were as¬ 
sumed to be equal to [pm]. The volume of an adult M 
balthica was estimated from the wet weight, assuming 
that 1 g wet weight corresponds to 1 cm^. The shape 
coefficient dm, used to relate volumetric shell length 
(SL; pm^) to volume (V; pm^) was (SL x 6^)^/ 
V=0.386. Finally, physiological rates were corrected 
for temperature, using an Arrhenius temperature of 
7596 K (M edulis'. Van Haren and Kooijman, 1993). 

3. Results 

3.1. Concentrations of bivalve larvae 

Larval concentrations in Lim^orden generally 
reached peak values in May with 117-241 larvae 
dm“^. At stations LF-LB and LF-SF, only single 
peaks were observed, while at station LF-NB, three 
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peaks were observed in May, June and August (Fig. 4). 
In the Sylt-Romo bight larval concentrations were 
slightly lower than in Limiiorden and the Delta area 
and consistently lower than in the western Wadden Sea. 
Due to the high sampling frequency, several smaller 
peaks were observed in the Sylt-Romo bight. However, 
a first large larval peak (80 larvae dm“ was observed 
at the beginning of May. A second peak (16 larvae 
dm“^) was observed in the first week of August and 
a third peak (13 larvae dm“^) at the beginning of 
September (Fig. 4). In the western Wadden Sea 
(WWS) the first clear peak (272 larvae dm“^) was 
noticed at the beginning of May, and the second peak 
at the end of June (392 larvae dm“^), the highest of all 
stations. Larvae were only occasionally found in winter 
and early spring samples. From the end of April to the 
end of October larval concentrations varied between 10 
and 100 larvae dm“^ (Fig. 4). In the Oosterschelde, 
peak values in the central (OS-C; 156 larvae dm“^) and 
western part (OS-W; 224 dm“^) were found in the third 
week of April, while peak values in the north (OS-N; 
103 dm“^) and eastern part (OS-E; 224 dm“^) were 
found in the third week of July. Concentrations of 
bivalve larvae in the Westerschelde (WS) were consis¬ 
tently lower than concentrations in the Oosterschelde. 
The maximum value was 16 larvae dm“^ (Fig. 4). 

To see whether larval concentrations depend on the 
water volume of the sampling areas, the total reproduc¬ 
tive output per basin was estimated from literature data 
on stock size, individual egg production and area sur¬ 
face (Tables 3, 4) and divided by the total water volume 
of each basin to obtain an egg concentration index 
(eggs dm“^). Next, the egg concentration index (eggs 
dm“^) was compared to the average larval concentra¬ 
tion (larvae dm“^) (Fig. 5). The measured larval con¬ 
centrations show a positive relation with the egg 
production index of the bivalve stock in each area 
(Tables 3, 4; Fig. 5). 

3.2. Food availability for bivalve larvae 

Among the 311 sampled phytoplankton taxa, 52 taxa 
were selected as suitable for larval ingestion (Table 5). 
The carbon content of this phytoplankton fraction var¬ 
ied from 2.0 x 10“’ to 2.1 x 10“' pg C cm“^ corres¬ 
ponding to an energy content of 9.4 x 10“^ to 9.5 x 
10“^ J cm“^. In general, possibilities for growth were 
maximal between May and September. At most sta¬ 
tions, the energy concentration of the phytoplankton 
was sufficient to meet maintenance costs (M; J larva“' 
d“^; Fig. 6). However, in the Sylt-Romo bight, larval 
maintenance costs were not reached for at least half of 


Table 3 

Main characteristics of study areas 


Area 

(km") 

(nh 

Volume 
(10‘" 1) 

Residence 
time (d) 

Limfjorden 

1575(‘> 

4.9(1) 

7.72 


L0gst0r Bredning (LF-LB) 




35<"> 

Skive Fjord (LF-SF) 




70(2) 

Nissum Bredning (LF-NB) 




15(2) 

Sylt-R0m0 bight (Sylt) 

400‘"> 

4.2 

1.69<'‘> 

19-29(5) 

Western Wadden 

1390<'^> 

3.3<"> 

4.59 

5-15(6) 

Sea (WWS) 





Delta area 





Oosterschelde 

350<*> 

7.7 

2.70<"> 


Western region 





(OS-W, OS-N) 





Central region (OS-C) 




140<"> 

Westerschelde (WS) 

300<‘“> 

8 

2.40 



References are indicated with numbers in brackets. 


References: (1) Dolmer and Frandsen, 2002; (2) Dolmer, 2002; (3) 
Strasser and Gunther, 2001; (4) Strasser, 2000; (5) Sprung et at, 2001; 
(6) Dame and Prins, 1998; (7) Ridderinkhof et al., 1990; (8) Bakker, 
1994; (9) Smaal et al, 2001; (10) Bouma et al, 2001. 

the year. Furthermore, the period in which a larva 
could grow was relatively short compared to other 
stations. At all stations, the most important contributor 
to total energy content of the phytoplankton was the 
category ‘other nanoplankton’ including dinoflagellates 
(Table 5) (94.4-99.5%) followed by diatoms (4.0- 
0.4%) and ciliates (0.0-1.5%). The average biovolume 
of the selected phytoplankton was 28 ± 20 pm^ (mean- 
± SE), ranging between 6.2 pm^ (station LF-LB) and 
92.7 pm^ (station WWS). In general, the potential 
assimilation rates (pa) ranged between >1% and 60% 
of the maximum assimilation rates (pAm) in the period 
between the beginning of April and the end of Sep¬ 
tember. Average values for this period were 18 ± 3 
(mean ± SE) to 26 ± 5% for the Limfjorden stations, 
7 ± 2% for the Sylt-Romo bight, 19 ± 2% for the west¬ 
ern Wadden Sea and 15 ± 3 to 25 ± 4% for the Oos¬ 
terschelde stations. 

In Limijorden, peak values of the phytoplankton en¬ 
ergy concentration occurred at the end of Jime and the 
beginning of July, varying between 7.0 x 10“^ J cm“’ 
(station LF-LB) and 9.5 x 10“^ J cm“^ (station LF-SF). 
According to the DEB-model, the larvae assimilated the 
energy content of the phytoplankton during these peaks at 
rates (pa) of 1.57 x 10“^ and 1.64 x 10“^ Jlarva“' d“ ^ 
respectively (Fig. 7). In the Sylt-Romo bight, phytoplank¬ 
ton energy concentration showed very distinct peak 
values at the beginning of May (8.4 x 10“^ J cm“^) 
and in June (6.4 x 10“^ J cm“^), while the energy con¬ 
centration was close to zero during the rest of the year. 
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Table 4 

Estimation of total potential bivalve reproductive output (eggs dm“^) per study area in 2000 



Species 

FW 

(10" kg) 

SL 

Ind. 

FW 

(g) 

Stock size 
(10'* ind.) 

Repr. 
output 
(lO" eggs 
female^') 

Total repr. 
output 
(lO*" 
eggs) 

Total egg 
density 
(eggs 
dm^^) 

Mean egg 
density* 
(eggs 
dm^^) 

Mean egg 
density* 
(eggs 
dm“^) 

Limfjorden 

Mytilus edulis 

400l" 

50<‘l 

10.8<"> 

37.02 

15.0l^> 

278 

3598 

23.7 

176305 

Sylt-Ri0m0 

Cerastoderma edule 

7.80<''l 

28 

7.0l"> 

1.12 

4.50l^> 

2.52 

149 

0.9 

6300 

Bight 

Crassostrea gigas^ 




0.00 F"! 

100<"> 

0.05 

3 

0.02 

125 


Macoma balthica 

3.64W 

16 

0.8<"> 

4.61 

0.45^1 

1.04 

61 

0.4 

2592 


Mytilus edulis 

1.48W 

50 

10.8(^> 

0.14 

15.0^1 

1.03 

61 

0.4 

2569 

Western 

Cerastoderma edule° 




21.3<"> 

4.50l^> 

48.0 

1047 

6.9 

34558 

Wadden 

Macoma balthica'' 




43.5l’> 

0.45l^> 

9.78 

213 

1.4 

7034 

Sea 

Mytilus edulis cult. 

39.0<*1 

50 

10.8 

3.61 

15.0® 

27.1 

590 

3.9 

19478 


Mytilus edulis subtid. 

33.7<^i 

50 

10.8 

3.12 

15.0® 

23.4 

510 

3.4 

16831 


Mytilus edulis intert 

2,15(10) 

50 

15.1 

0.15 

15.0® 

1.16 

25 

0.2 

832 

Delta area 











Oosterschelde 

Cerastoderma edule 


28(11) 

7.0<"> 

7.61<’> 

4.50® 

17.1 

634 

4.2 

48921 


Crassostrea gigas" 

150O^> 


150 

1.00 

100<"> 

50.0 

1852 

12.2 

142857 


Crassostrea gigas 

6.500^1 


150 

0.043 

100<"> 

2.17 

80 

0.5 

6190 


Macoma balthica 




2.04^1 

0.45® 

0.46 

17 

0.1 

1312 


Mytilus edulis^ 




3.74<‘'> 

15.0® 

28.1 

1039 

6.8 

80143 

Westerschelde 

Cerastoderma edule 




1.02^1 

4.50l^> 

2.29 

109 

0.6 

7643 


Macoma balthica 




4.42<'*> 

0.450<^1 

0.995 

47 

0.3 

3318 


Only the most important bivalve species are included. Estimates of the reproductive output are based on stock estimates for 2000, except for the 
Sylt-Romo bight where estimates were based on surveys in the mid-1990s. When direct stock size estimates were not available, stock size was 
estimated assuming a species-specific standard shell length. It is assumed that 50% of the bivalve stock consists of females. FW= fresh weight, 
SL=shell length, ind.=individual, repr.=reproductive, cult.=subtidal culture plots, subt. = subtidal, intertid.=intertidal. References are indicated 
with numbers in brackets. 

References: (1) Kristensen and Hoffinann, 2004; (2) Williams and Ens, unpubl. results; (3) Honkoop et al., 1999; (4) Gatje and Reise, 1998; (5) 
Zwarts and Blomert, 1992; (6) Reise, 1998; (7) Kamermans et al, 2004; (8) Smaal et al, 2001; (9) Kater and Kesteloo, 2003; (10) Kater and Den 
Os, 2001; (11) Dekker and De Bruin, 2001; (12) Perdon and Smaal, 2000. 

“ Excluding hard substrates. 

Assuming that the stock in western Wadden Sea is 50% of the stock in the Dutch Wadden Sea (Brinkman, pers. comm.). 

* Averaged over the length of the spawning season (1 April-30 Sept, 152 days) for comparison with larval densities. 


Assimilation rates (pa) for these peak values were 
1.7 X 10“^ and 2.7 x 10“^ J larva"' d"', respectively 
(Fig. 7). In the western Wadden Sea, the energy concen¬ 
tration fluctuated over the year, with peak values in 
February (7.8x10“^ J cm"^), April (6.1x10“^ J 
cm“^) and June/July (4.9 x 10“^ J cm“^), yielding as¬ 
similation rates Pa of 1.55x10“^, 1.46x10“^ and 
2.5 X 10“^ J larva"'d"', respectively (Fig. 6). 

In the Delta area, peak values were mainly found in 
May at stations OS-N and OS-E (~2.8 x 10"^ J cm"^), 
at end of Jime at all three stations (4.0^.6x 10"^ J 
cm"^) and in August at stations OS-W and OS-N (2.1- 
4.2 X 10"^ J cm"^). Assimilation rates of larvae were 
1.33-1.34X 10"^ J larva"' d"' (May, OS-N and OS- 
E), 3.01-3.97 J larva"' d"' (June, all three stations), 
and 3.90-2.87x 10"^ J larva"' d"' (August, OS-W 
and OS-N). Due to the absence of phytoplankton data 
for stations OS-C and WWS, no assimilation rates 
could be estimated. 


4. Discussion 

4.1. Assumptions and model limitations 

The results of the present study suggest that larvae in 
north-west European waters are food limited through¬ 
out the year. It has to be stressed, however, that all 
estimates are based on assumptions. Estimates of larval 
abundance were made imder the assumption that larvae 
are distributed uniformly over the water column, so any 
vertical movements of the larvae (Verwey, 1966) were 
ignored. Average larval concentrations in the water 
column may, hence, have differed from those at the 
surface, where the sampling took place. In addition, it 
was assumed that larvae would not be size-selective and 
would ingest all phytoplankton and heterotrophic cili- 
ates <20 pm, except for species possessing spines or 
forming chains. However, some species of larvae are 
known to feed selectively. For example, larvae of the 
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Fig, 5. An estimated index of average bivalve egg concentrations 
(eggs dm“^) plotted against the average number of larvae measured in 
each area (larvae dm“^; mean ± SE). Egg concentrations are based on 
total yearly egg production (Table 4) and the water volume of the 
sampled areas (Table 3). Results of De Vooys (1999) (station WWS in 
earlier years) have been included for comparison (De Vooys’ data on 
Mytilus edulis and Cerastoderma edule biomass were converted to 
egg numbers according to Table 4). 

oyster Crassostrea virginica strongly select for 2—A |Tm 
food particles in food suspensions, and for much larger 
(22-30 pm) food material in the presence of a large 
(>10 pm)-cell dinoflagellate bloom (Baldwin, 1995), 
while Mytilus edulis veligers had the highest clearance 
rates at particle sizes of 2.5 to 3.5 pm and did not ingest 
particles <1 and >8-9 pm (RiisgM and Randlov, 
1980). The volumes of the prey cells were estimated 
from the literature, and their specific size was assumed 
to be constant. The nutritional value (J pm“ ^) of phy¬ 
toplankton was also assumed to be constant, although 
this can differ between and within species (Walne, 
1963; Baldwin and Newell, 1995). 

It was assumed that the Dynamic Energy Budget 
(DEB) model for Macoma balthica larvae (Philippart et 
al., ms. in prep.) would provide good estimates of the 
energetic requirements. This probably holds for main¬ 
tenance costs (M) (Kooijman, 2000; Van der Veer et al., 
2001), but not for ingestion rates. Van der Veer et al. 
(2001) showed that various flatfish species differ espe¬ 
cially with respect to their surface area-specific maxi¬ 
mum ingestion rates {Jxm}- For larvae of various 
bivalve species this may mean that ingestion rates differ 
and, hence, that the degree of food limitation for other 
species than M balthica may have been under- or 
overestimated. In addition, the energy budget calcula¬ 
tions were made for a 200-pm larva, although in the 
field larvae may range from 50 to 400 pm (Pulfrich, 
1997; pers. obs.). According to the DEB-model, a 
smaller larva (e.g., 150 pm) and a larger one (e.g., 
250 pm) require 56% and 156%, respectively, of the 


energy required by a 200-pm larva (Fig. 7), which 
corresponds to the observation that smaller organisms 
generally have a higher demand per unit of body mass 
(e.g. Banse and Mosher, 1980; Banse, 1982). Finally, in 
the DEB-model itself, many assumptions have been 
made that are discussed in detail by Philippart et al. 
(ms. in prep.). For example, the assimilation efficiency 
in de DEB-model is constant (0.56%; Van Haren and 
Kooijman, 1993), but may in reality show great varia¬ 
tion (33 to 92% in mixed zooplankton; Conover, 1966). 
Altogether, numerous xmcertainties exist in the esti¬ 
mates of the parameters. Confidence intervals could 
be provided for the results, based on the imcertainties 
of the parameters, but that would only suggest a cer¬ 
tainty of the data that is not realistic. 

4.2. Larval concentrations 

In the current study, the observed average larval 
concentrations (larvae dm“^) were positively related 
to the egg concentration index (eggs dm“^) (Fig. 5), 
i.e. larval concentrations appeared to depend on the size 
of the adult stocks and the water volume of the estuary. 
To keep it simple, this analysis did not include the water 
exchange of the basins with the open sea, resulting in 
offshore transport of larvae (Roegner, 2000). The egg 
concentration index can probably be improved if more 
information on the various stock sizes was available. In 
the western Wadden Sea, for example, the deeply bur¬ 
ied Mya arenaria was not taken into accoimt because 
only limited information exists on its distribution, al¬ 
though in parts of the western Wadden Sea it may 
constitute 7-70% of the total bivalve biomass (Dekker 
and De Bruin, 2001). Similarly, Ensis americanus in 
the Sylt-Romo bight was not taken into account be¬ 
cause only limited information was available on its 
abundance and distribution (Armonies and Reise, 
1999). In addition, the population size of Crassostrea 
gigas in the Sylt-Romo bight may have been under¬ 
estimated, since it increased from 1 million individuals 
in 1995 (Reise, 1998) to 60 million individuals in 2003 
(Diederich et al, 2005). 

The data of De Vooys (1999), who measured larval 
concentrations at station WWS in earlier years, confirm 
the positive relation between the egg production index 
and the measured larval abundance at low egg concen¬ 
trations, i.e. at small adult stock (Fig. 5). The egg 
production index was estimated by converting estimates 
on the fresh weight of the mussel and cockle stocks (De 
Vooys, 1999) into egg numbers according to Table 4. At 
extremely high bivalve stocks (as occurred around 1981 
due to simultaneous successful recmitment in all im- 
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Table 5 

Individual cell volumes (nm^) of phytoplankton taxa that can theoretically be ingested by bivalve larvae on basis of their size and morphology (see 
methods section), sampled in Limfjorden, Sylt-Romo bight, western Wadden Sea and Delta area in 2000 


Location 

0”“- 

Species 

s-"* 

Calculation of 

biovolume 

V(pm2) 

pgC cell ' 

Reference 

Limfjorden 

0 

Autotrophic picoplankton indet. 

6 

ir/6 X (1)^ 

0.52 

0.25 

a 

Limfjorden 

D 

Centric diatom <10 pm 

1 

17/4 X (5)^ X 2.5 

49.09 

7.24 

b 

Sylt-Romo bight 

D 

Centric diatom 0-2 

1 

17/4 X (1)^ X 0.5 

0.39 

0.19 

b 

Sylt-Ramo bight 

D 

Centric diatom 0-3 

1 

17/4x(1.5)^x0.75 

1.33 

0.47 

b 

Sylt-Rama bight 

D 

Centric diatom 0-5 

1 

17/4 X (2.5)^ X 1.25 

6.14 

1.50 

b 

Lim^orden 

D 

Centric diatom 10-20 pm 

1 

17/4 X (15)^ X 7.5 

1325.36 

88.05 

b 

Sylt-Rama bight 

D 

Centric diatom 1-3 

1 

17/4 X {if X 1 

3.14 

0.90 

b 

Sylt-Rama bight 

D 

Centric diatom 2 

1 

17/4 X {if X 1 

3.14 

0.90 

b 

Sylt-Rama bight 

D 

Centric diatom 2x4 

1 

17/4 X {if X 1 

3.14 

0.90 

b 

Sylt-Rama bight 

D 

Centric diatom 5 

1 

17/4 X (5)^ X 2.5 

49.09 

7.24 

b 

Sylt-Rama bight 

D 

Centric diatom 5-10 

1 

17/4 X (7.5)^ X 3.75 

165.67 

18.20 

b 

Lim^orden 

0 

Choanoflagellates <5 pm 

5 

17/6 X (1.25)^ X 2.5 

2.05 

0.81 

b 

Limfjorden 

0 

Choanoflagellates 5-10 pm 

5 

17/6 X Q.lSf X 7.5 

55.22 

13.82 

b 

Limfjorden 

0 

Chrysochromulina spp. <5 pm 

6 

17/6 X (2.5)^ 

8.18 

2.66 

b 

Lim^orden 

0 

Chrvsochromulina spp. 5-10 pm 

6 

17/6 X (7.5)^ 

220.89 

45.71 

b 

Lim^orden 

C 

Ciliates spp. 0-10 pm 

5 

17/6 X (2.5)^ X 5 

16.36 

2.70 

a 

Limfjorden 

C 

Ciliates spp. 11-20 pm 

5 

17/6 X (7.5)^ X 15 

441.79 

72.90 

a 

Limfjorden 

0 

Coccoid cyanobacteria spp. 

5 

17/6 X (1)2 X 2 

1.05 

0.45 

a 

Delta area 

0 

Cryptomonadales 

6 

17/6 X (4)2 

33.51 

8.98 

b 

Delta area 

0 

Cryptophyceae sp. 

5 

17/6 X (1.25)2 X 2.5 

2.05 

0.81 

b 

Limfjorden 

0 

Cryptophyceae spp. <5 pm 

5 

17/6 X (1.25)2x2.5 

2.05 

0.81 

b 

Limfjorden 

0 

Cryptophyceae spp. 10-15 pm 

5 

17/6 X (6.25)2x12.5 

255.66 

51.86 

b 

Lim^orden 

0 

Cryptophyceae spp. 15-20 pm 

5 

17/6 X (8.75)2x17.5 

701.54 

123.92 

b 

Lim^orden 

0 

Cryptophyceae spp. 5-10 pm 

5 

17/6 X (3.75)2 ^ .7 5 

55.22 

13.82 

b 

Western Wadden Sea 

D 

Cymatosira belgica 

- 

- 

300 

59.53 

c 

Delta area 

D 

Delphineis minutissima 

- 

- 

27 

7.45 

d 

Western Wadden Sea 

0 

Dinoflagellate 

5 

7t/6x(7)2x 10 

256.56 

52.02 

e 

Sylt-Rama bight 

0 

Dinoflagellate sp. 

6 

7t/6x(7)2x 10 

256.56 

52.02 

e 

Sylt-Rama bight 

0 

Dinoflagellate sp. >3 pm 

5 

77/6 X (1.5)2 X 3 

3.53 

1.29 

b 

Sylt-Rama bight 

0 

Dinoflagellate sp. >4 pm 

5 

77/6 X (2)2 X 4 

8.38 

2.71 

b 

Sylt-Rama bight 

0 

Dinoflagellate sp. 0-3 pm 

5 

77/6 X (0.75)2 X 1.5 

0.44 

0.21 

b 

Lim^orden 

0 

Dinoflagellate sp. 10-20 pm 

5 

77/6 X (7.5)2 ^ J5 

441.79 

72.90 

b 

Sylt-Rama bight 

0 

Dinoflagellate sp. 5 pm 

5 

77/6 X (2.5)2 X 5 

16.36 

4.84 

b 

Sylt-Rama bight 

0 

Dinoflagellate sp. 7 pm 

5 

77/6 X (3.5)2 ^ y 

44.90 

11.56 

b 

Lim^orden 

0 

Ebria sp. (10-20 pm) 

3 

17/12 X (8)2 X (4+7) 

181.34 

38.55 

b 

Western Wadden Sea 

D 

Eunotogramma dubium 

- 

- 

539 

44.52 

f 

Sylt-Rama bight 

D 

Eunotogramma dubium 

- 

- 

539 

44.52 

f 

Western Wadden Sea 

0 

Flagellates <10 pm 

- 

- 

60 

14.84 

c 

Limfjorden 

0 

Flagellates <5 pm 

5 

17/6 X (1.25)2x2.5 

2.05 

0.81 

b 

Lim^orden 

0 

Dinoflagellate sp. < 10 pm 

5 

17/6 X (2.5)2 X 5 

16.36 

4.84 

b 

Western Wadden Sea 

0 

Flagellates >10 pm 

- 

- 

612 

110.15 

c 

Limfjorden 

0 

Flagellates 5-10 pm 

5 

17/6 X (3.75)2 ^ .7 5 

55.22 

13.82 

b 

Delta area 

0 

Heterocapsa minima 

- 

- 

125 

27.97 

d 

Western Wadden Sea 

0 

Heterocapsa rotundata 
(syn. Katodinium rotundatum) 

3 

17/12 X (6)2 X (6.5+6) 

117.81 

26.57 


Sylt-Rama bight 

0 

Heterocapsa sp. 

5 

17/6 X (5.25)2 X 10.5 

151.53 

33.02 

b 

Delta area 

0 

Heterocapsa triquetra 

- 

- 

1953 

299.83 

d 

Lim^orden 

0 

Heterotrophic flagellates 

- 

- 

60 

14.84 

c 

Sylt-Rama bight 

D 

Katodinium sp. 

- 

- 

151.53 

33.02 

b 

Delta area 

0 

Leucocryptos marina 

3 

77/12 X (7.5)2 ^ ^.7 5^ 

316.61 

62.37 

a 

Lim^orden 

0 

Leucocryptos sp. (10-15 pm) 

3 

17/12 X (6)2 X (10+6) 

150.80 

32.88 

b 

Lim^orden 

0 

Leucocryptos sp. (5-10 pm) 

3 

17/12 X (3)2 X (5+3) 

75.40 

18.08 

b 

Lim^orden 

C 

Lohmaniella oviformis 

5 

77/6 X (15.5)2 X 16.5 

2075.61 

342.48 

a 

Limfjorden 

0 

Micromonas pusilla 

5 

77/6 X (1)2x2 

1.05 

0.45 

b 

Delta area 

0 

Nephroselmis sp. 

4 

17/4 X 6.5 X 4 X 2 

40.84 

10.65 

a 

Western Wadden Sea 

D 

Pennate diatoms indet. 

- 

- 

895 

65.39 

c 






O.G. Bos et al. /Journal of Sea Research 55 (2006) 191-206 


201 


Table 5 {continued) 


Location 

o„.p 

Species 


Calculation of 

biovolume 

V(pm2) 

pgC cell ‘ 

Reference 

Lim^orden 

D 

Pennate diatoms indet. 

- 

- 

895 

65.39 

c 

Delta area 

0 

Peridiniaceae 

- 

- 

4096 

568.15 

d 

Western Wadden Sea 

0 

Phaeocystis flagellate 

- 

- 

65 

15.91 

c 

Sylt-Romo bight 

0 

Phaeocystis globosa 

- 

- 

65 

15.91 

c 

Delta area 

0 

Phaeocystis sp. 

- 

- 

65 

15.91 

c 

Delta area 

0 

Prasinocladus marinus 

5 

7t/6 X (6)^ X 10 

188.50 

39.87 

b 

Limfjorden 

0 

Prasinophyceae spp. 

5 

7i/6 X (6)^ X 10 

188.50 

39.87 

b 

Sylt-Romo bight 

0 

Prasinophyceae spp. 

5 

7t/6x(6)2x 10 

188.50 

39.87 

b 

Lim^orden 

0 

Prorocentrum balticum 

- 

- 

375 

72.18 

g 

Lim^orden 

0 

Prorocentrum minimum 

5 

77/6 X (12.5)^ X 18 

1472.62 

235.00 

a 

Limfjorden 

0 

Prorocentrum triestinum 

5 

ir/6x(8)^x 10 

335.10 

65.50 

a 

Delta area 

0 

Pyramimonas sp. 

2 

ir/12x(10)^x 15 

392.70 

75.11 

b 

Western Wadden Sea 

0 

Rhodomonas spp. 

- 

- 

221 

45.73 

c 

Delta area 

0 

Small algae indet. 

6 

77/6 X (2.1)^ 

4.85 

1.69 

d 

Limfjorden 

c 

Strombidium sp. (<10 pm) 

2 

77/12 X (3)^ X 5 

11.78 

1.94 

h 

Limfjorden 

c 

Strombidium sp. (10-20 pm) 

2 

77/12 X (9)2 X 15 

318.09 

52.48 

h 


D=diatoms, C=ciliates, 0=other nanoplankton (including dinoflagellates). Carbon content of diatoms is calculated according to Strathmann 
(1967), of ‘other nanoplankton’, according to Verity (1992) and of ciliates according to Putt and Stoecker (1989) (see ‘Materials and methods’). 
References: (a) Kuylenstiema and Karlson, 1998; (b) Tomas, 1997; (c) Philippart et al., 2000; (d) Koeman, pers. comm., 2003; (e) Van Iperen, 
pers.comm., 2003; (f) Philippart pers. comm., 2003; (g) Plymouth Mar. Lab.; (h) Ichise et al., 2004. 

’* Geometric shapes and biovolume calculations as described in Hillebrand et al, 1999: (1) cylinder; (2) cone; (3) cone + half sphere; (4) elliptic 
prism; (5) prolete spheroid; (6) sphere. 


portant bivalve species in 1979; compare Beukema et 
al., 2001) however, De Vooys (1999) measured low 
larval concentrations, probably due to the important 
filtration activity of the adults. Furthermore, De 
Vooys (1999) showed a clear year-to-year variation in 
larval concentrations in a single area (western Wadden 
Sea), which demonstrates that geographical compari¬ 
sons of larval concentrations based on different studies 
from different years must be treated with great caution. 

Maximum larval concentrations in the current study 
varied from 17 (station WS) to 392 larvae dm“^ (station 
WWS), with intermediate maximum values of 100-220 
larvae dm“^ at the other stations. The low larval con¬ 
centrations in the Westerschelde may be attributed to a 
combination of relatively low bivalve stocks and a 
relatively large water volume (Fig. 5). In the Sylt- 
Romo bight, at similarly low bivalve stocks (Table 4), 
larval concentrations were higher, because the eggs 
were probably less diluted (Fig. 5). Larval numbers 
were probably underestimated in the Sylt-Romo bight 
because the mesh size of the sieve was larger than in the 
other study areas (Table 1). In addition, the adult stock 
may have been imderestimated. In the western Wadden 
Sea, the highest larval concentrations were observed, 
which was attributed to a relatively low water volume 
and high bivalve stocks (Fig. 5). In the Oosterschelde 
and Limfjorden, larval concentrations were of the same 
magnitude. In the Oosterschelde, bivalve stocks were 
larger, but the dilution of eggs was also greater. Since, 


in general, the average benthic biomass is limited by the 
primary productivity of the system (Flerman et al., 
1999), it may be that larval concentrations are indirectly 
determined by the primary production, via adult bio¬ 
mass and adult reproductive output. 

The study areas differed not only in the intensity of 
the larval peaks, but also in their number. The succes¬ 
sive larval peaks may have been the result of the 
different timing of spawning of bivalve species. Single 
larval peaks were observed in Limfjorden, while in the 
other areas multiple larval peaks were observed. This 
may be an effect of the sampling frequency (Table 1), 
but also of differences in physical and hydrodynamic 
characteristics of the areas (Table 3). Limfjorden is a 
subtidal area, whereas tidal flats dominate the other 
areas. In general, the water temperature increase in 
subtidal areas is spatially more uniform than in inter¬ 
tidal areas (Raffaelli and Hawkins, 1996), leading to a 
more synchronised spawning in subtidal areas. Further¬ 
more, Limijorden is dominated by a single bivalve 
species, M. eduUs (Dolmer et al., 1999), while in the 
other areas 3-5 bivalve species are abimdant (Table 4). 
Therefore the single larval peaks observed at stations 
LF-LB and LF-SF may have resulted from a synchro¬ 
nised spawning event of M. edulis, while the multiple 
peaks measured in the other areas (Fig. 4) may have 
resulted from multiple spawning by multiple species. 
Finally, the duration of peaks may vary between areas 
due to effects of water temperature (Drent, 2002) and 
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Fig, 6. Assimilation rate of a bivalve larva (p*; J larva^' d“'; middle line) based on phytoplankton energy content and larval energy needs. The 
upper line shows maximum assimilation rate (pAm; J larva^' d^*) leading to maximum larval growth, and lower line shows maintenance costs (M; J 
larva^' d *), ut which growth is not possible. Rates are calculated with a DEB-model for a 200 pm Macoma balthica larva and are corrected for 
temperature. 


nutrition (Pechenik et al., 2002) on larval development 
rate. 

4.3. Food limitation 

The above results show that in the period between 
April and September 2000, larvae were able to achieve 


an average assimilation rate Pa between 7% and 25% of 
the maximum assimilation rate pAm, with a maximum 
of 60%. This suggests that a representative larva, with a 
length of 200 pm, is food-limited for most of its life, 
irrespective of the study site. Such low assimilation 
rates may greatly increase the time a larva needs to 
complete development into a competent larva, with the 



Fig. 7. DEB-model predictions of energy intake (J larva^' d^*) as a function of food concentration (functional response; cells Isochrysis 
galbana cm“^) for different larval sizes (150, 200 and 250 pm). From top to bottom, lines indicate maximum ingestion rate (Jxm)i ingestion 
rate (Jx). maximum assimilation rate (pAm)i assimilation rate (pa), and maintenance costs (M). The vertical line indicates the half-saturation 
constant Xi^. 
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risk of being transported offshore (Young et al., 1998; 
Roegner, 2000) or of predation by e.g. benthic filter 
feeders (Cowden et al., 1984; Andre and Rosenberg, 
1991). However, by definition, larvae will never attain 
maximum growth rates and will therefore always be 
food-limited to a certain extent, since the pAm is the 
upper limit of an asymptotic function (Fig. 3). A 
relevant question is therefore at what food levels a 
larva is still able to grow and develop successfully. In 
a laboratory experiment, at 15 °C, M. balthica larvae 
were able to develop successfully to metamorphosis at 
food concentrations of 4.0 x 10^ cells Isochrysis gal- 
bana cm“^ (Bos et al., subm. ms.). For a 200-pm 
larva, this would correspond to a Pa of 29% of the 
PAm- Furthermore, growth rates were high and devel¬ 
opment was fast at food concentrations of 80 x 10^ 
cells cm“^ (Bos et al., subm. ms), corresponding to a 
Pa of 89% of PAm for a 200-|am larva at 15 °C. In 
another experiment, M. balthica larvae reared in nat¬ 
ural seawater were mostly not able to complete meta¬ 
morphosis, while control larvae fed with 50 x 10^ I. 
galbana cells cm“^ always successfully completed 
metamorphosis (Bos et al., subm. ms). Together, 
these experiments suggest that the observed food 
levels in the field were frequently too low for bivalve 
larvae to complete metamorphosis. 

Other studies in coastal waters that aimed at detect¬ 
ing food limitation of larvae have yielded contrasting 
results. In a Danish embayment, mussel M edulis 
larvae reared in natural seawater were generally 
food-limited, compared to those reared in emiched 
natural seawater (Fotel et al., 1999). Jorgensen 
(1981), however, reported growth rates of M. edulis 
veligers estimated irom field samples similar to growth 
rates under optimal conditions in the laboratory. Thus 
the literature suggests food limitation in at least some 
coastal areas. Limited food availability in our study 
areas is probably not caused by competition between 
larvae (Strathmann, 1996; Bos et al., ms. in prep.), but 
by competition of larvae and the benthic bivalve stock 
(Cowden et al., 1984; Andre and Rosenberg, 1991), a 
match/mismatch situation of larvae and their food 
source (Philippart et al., 2003; Bos et al., subm. ms.) 
as has been reported for fish larvae (Cushing, 1990), 
or by competition with other planktonic species (Han¬ 
sen et al., 1994). 

In conclusion, the results suggest that larvae in 
north-western European waters are food limited 
throughout the year under the assumptions made in 
the present study. Although the assimilation rates at 
most study sites appear to be high enough to cover 
the maintenance costs of the larvae, it is not clear 


whether they can support the metamorphosis of veliger 
larvae into settlers. As a consequence, the link between 
the degree of larval food limitation and bivalve recruit¬ 
ment remains to be elucidated (Olson and Olson, 1989). 
Further research should focus on the determination of 
natural larval food sources, their nutritional value, and 
their effect on larval growth, development and survival 
rates. Also factors that influence larval food availability, 
such as competition for food with other planktonic 
species or with the adult stock, and the behaviour of 
larvae in response to patches of food (Metaxas and 
Young, 1998) should be clarified. 
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